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Cognitive Radio with Partial Channel State 
Information at the Transmitter 

Abstract 

In this paper, we present the cognitive radio system design with partial channel state information 
known at the transmitter (CSIT). We replace the dirty paper coding (DPC) used in the cognitive radio with 
full CSIT by the linear assignment Gel'fand-Pinsker coding (LA-GPC), which can utilize the limited 
knowledge of the channel more efficiently. Based on the achievable rate derived from the LA-GPC, 
two optimization problems under the fast and slow fading channels are formulated. We derive semi- 
analytical solutions to find the relaying ratios and preceding coefficients. The critical observation is that 
the complex rate functions in these problems are closely related to ratios of quadratic form. Simulation 
results show that the proposed semi-analytical solutions perform close to the optimal solutions found by 
brute-force search, and outperform the systems based on naive DPC. Asymptotic analysis also shows 
that these solutions converge to the optimal ones solved with full CSIT when the /T-factor of Rician 
channel approaches infinity. Moreover, a new coding scheme is proposed to implement the LA-GPC in 
practice. Simulation results show that the proposed practical coding scheme can efficiently reach the 
theoretical rate performance. 

I. INTRODUCTION 

As the demand for high data rate steadily increases, efficient spectrum usage becomes a 
critical issue. Recent measurements from the Federal Communications Commission (FCC) have 
indicated that 90 percent of the time, many licensed frequency bands remain unused. A radio 
technology that promises to solve such problems is the cognitive radio (CR) [1]. This technology 
is capable of dynamically sensing and locating unused spectrum segments in a target spectrum 
pool, and communicating via the unused spectrum segments without causing harmful interference 
to the primary users. The primary users are defined as those of existing commercial standards. 
If a primary user demands the channel, the CR user should vacate and find an alternative one. 
Assuming that there is perfect channel state information at the transmitter (CSIT), Devroye et. 
al. in [2] proposed the concept of interference mitigation CR in which the CR and primary 
users can simultaneously transmit at the same time and frequency bands. The CR transmitter 
not only transmits its own signal but also relays the primary users'. This is the key to keeping 
the primary users' transmission rates unchanged under the coexistence of CR users since the CR 
users' signals may interfere with the primary users'. For the CR receiver, signals from the primary 
transmitters and those relayed by the CR transmitter are interferences. These interferences are 
non-causally known at the CR transmitter and can be precoded with dirty paper coding (DPC) 
[3] such that the CR receiver can decode as if these interferences do not exist. The validity of 
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the non-causally known interferences has been discussed in [4]. DPC is a promising precoding 
technique for cancelling arbitrary interferences perfectly known only at the transmitter but not at 
the receiver. Devroye et. al. [2] then derived the achievable rate of the DPC based interference 
mitigation CR system. In [4], Jovicic and Viswanath further proved that the DPC based CR 
system indeed achieves the capacity. 

Due to limited feedback bandwidth and the delay caused by channel estimation, it is more 
practical to consider the case with only partial CSIT. Unfortunately, the interference-free rate 
achieved by DPC heavily relies on perfectly known CSIT. It has been observed that naively 
applying DPC in fading channels with partial CSIT may cause a significant performance loss 
[5]. In this paper, we focus on the CR system design for the cases where only the fading 
channel statistics are available at both the primary and CR transmitters. The main contributions 
of this paper are as following. We use a coding scheme more general than the DPC, termed 
the linear-assignment Gel'fand-Pinsker coding (LA-GPC) [6] to derive the achievable rate of 
the CR channel. Based on the achievable rate, we then propose design methods to optimize the 
performance of the CR system with only the statistics of CSIT for both fast and slow fading 
scenarios [7]. The goal is to keep the primary user's rate unchanged while maximizing the 
CR user's rate. The considered problems are non-convex and are almost analytically intractable. 
Therefore, we provide semi-analytical solutions to the optimal relaying ratio of the CR transmitter 
in fast and slow fading channels, respectively. These solutions are obtained by expressing the 
primary user's rate constraint in the ratio of quadratic form, then applying the techniques of 
moments of quadratic forms and the deterministic approximations of probability inequalities [8]. 
To optimize the CR user's own rate, we also provide semi-analytical solutions to the precoding 
coefficients in fast and slow fading channels, respectively. In slow fading channels, we show that 
the complex objective function can be recast as a compound function of the incomplete gamma 
function from which the precoding coefficient can be easily numerically solved. Asymptotic 
analysis is also given to show that these solutions converge to the optimal ones solved with 
full CSIT when the ^-factor of Rician channel approaches infinity. This result verifies the 
correctness of the proposed methods. Simulation results show that the proposed semi-analytical 
methods perform well compared to the optimal solutions found by brute-force search under 
various channel conditions. In addition to the theoretical results based on unstructured random 
Gaussian codebooks, we also use nested-lattice codebooks with lattice decoding to implement 
the precoding in practice [9]. Simulation results show that the nested-lattice coding scheme can 
efficiently approach the theoretical achievable rate. 

This paper is organized as follows. Section |II] gives the system model of the CR system with 
partial CSIT. Review of the LA-GPC and its application to the partial CSIT CR system are 
also given in this section. The proposed parameter design for fast and slow fading channels 
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are discussed in Section IIII-AI and Section IIII-Bl respectively. The asymptotic analysis of the 
parameters is given in Section |IVl A nested-lattice based LA-GPC scheme using the proposed 
design parameters is illustrated in Section |Vl Simulation results are given in Section |VIl Finally, 
Section IVIII concludes this paper. 

II. Background and System Model 

In this section we first introduce the system model in consideration. We then review the 
problems of DPC and LA-GPC in fading channels. After that we will derive the rate formula 
of the fading cognitive channel with LA-GPC Q 

A. The fading cognitive channel 

The cognitive channel coined by [2] is shown in Fig. [TJ The channel gains between pri- 
mary transmitter and receiver, primary transmitter and CR receiver, CR transmitter and primary 
receiver, and CR transmitter and receiver are denoted by //ii,//2i,^^i2, and H22, respectively. 
The primary and CR transmitted signals are Xp and Xc with transmit power constraint Pp and 
Pc, respectively. The unidirectional arrow from the primary transmitter to the CR transmitter 
means that the CR transmitter unilaterally knows the primary user's codewords non-causally. The 
feasibility of this assumption is discussed in [4]. The signals at the primary and CR receivers 
are 



Hn Hn 
H21 H22 



1 





1 



Xr 



+ 



(1) 



where Xc is the CR user's own signal after being precoded. The CR transmitter relays the 
primary user's signal to maintain the performance of the primary link. The relaying ratio ai 
is the percentage of CR's transmit power used for relaying the primary signal. Thus the factor 
a/ aiPc/Pp is to ensure that the relaying power is aiPc- The circularly symmetric complex 
additive white Gaussian noises (AWGN) at the primary and CR receivers are denoted by Zp ~ 
C3V;(0,o|^) and ~ C5V:(0,o| ), respectively. 

In the following discussion we assume that the primary receiver knows Hn perfectly, but 
only the statistics of Hn are known at primary transmitter. Similarly, the CR receiver knows all 
four channel gains perfectly, but only the statistics of them are known at the CR transmitter. We 



*In this paper, the superscript (.)^ and (.)^ denote the transpose and complex conjugate, respectively. Letters without and 
with underlines denote variables and vectors, respectively. If a letter is Italic capital, the corresponding variable or vector is 
random. Boldface capital letters denote deterministic matrices. A(!, j) denotes the entry at the rth row and the jth column of the 
matrix A. The covariance matrix of the random vector [AS]^ is denoted by Sa.b- |A| and \a\ represent the determinant of the 
square matrix A and the absolute value of the scalar variable a, respectively. l„ is the identity matrix with dimension n. I{A;B) 
denotes the mutual information between A and B. The trace of A is denoted by tr(A). The diagonal matrix whose diagonal 
entries are a is denoted by diag(a). All the logarithm operations are of base 2 such that the unit of rates is in bit. 
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assume that the primary and CR users are subject to slow or fast fading channels simultaneously. 
In this paper the four channels Hij, 1 < z, j < 2, are assumed to be Rice distributed with mean and 
variance denoted by /Jij and o?-, respectively. We also let \/Jij\'^ + ojj = 1 to keep the calculation 
of the received signal-to-noise ratio (SNR) simple. 

B. LA-GPC in fading channels 

Since only the statistics of the fading channels are known at the transmitter, the CR user 
cannot use the classical DPC to perfectly cancel the interference due to Xp. Instead, the DPC is 
replaced by the LA-GPC which is more suitable for this task. For brevity, the CR user's channel 
is recast as 

Y =H^X + HsS + Z, (2) 

with 



Hs = H21 + y ^^22, = H22, Y = Ys,X = Xc,S = Xp, and Z = Z„ (3) 

where X, S, and Z are assumed to be independent complex Gaussian random variables with zero 
mean and variances (powers) P^, P^, and = "^Zj' respectively. The fading channel gains Hj^ 
and Hs of the signal and interference, respectively, are known perfectly at the receiver, but only 
their statistics are known at transmitter. For illustration purpose, we first focus on the fast fading 
channel in this subsection. The capacity of this channel can be modified from [6] as 

sup E{R{H,,H,)}, (4) 
f{u\s)J{-) 

where 

R{H^, Hs) = I{U;Y\H^,Hs) - I{U;S) , (5) 

U is an auxiliary random variable with distribution obtained via the conditional distribution 
/ (mI*), and /(■) is a deterministic function such that X = /(f/, S). We select f{U, S) — U — 
which makes, 

U=X + a2S. (6) 

Such strategy function selection is called the linear assignment strategy and a2 is the preceding 
coefficient. With perfect CSIT, we can choose a2 as the linear minimum mean square error 
(MMSE) filter coefficient ttc = \Hjc\^Px/ {\Hx\^Pjc + Pz) to estimate X from the interference-free 
channel Y = H^X + Z, and the interference-free rate is achievable. Costa named the LA-GPC 
with a2 = CXc as DPC [3]. However, in our setting, the selection of the filter a2 must rely only 
on the channel statistics. Thus the MMSE filter used in Costa's DPC cannot be used here. 
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C. Achievable rate of CR with LA-GPC 

Before further discussing the primary and CR users' performances in fast and slow fading 
channels, the rate functions with given channel gains should be derived first. Assuming that 
the linear-assignment strategy is used, we have U = + CHiXp from ([3]) and The mutual 
information with given channel realizations can be computed as [10] 



R{h22,h2i)=I{U;Y,,\H22 = h22.H2i=h2i)-I{U;S) 

( 



l0g\^i\h22\W,^ + \h2l + ^^h22\^Pp + ol)al ) -l0g(|Et;,Kj), (7) 



where oj^ = {I —ai)Pc, and ^u,Ys is 

h22<sl+a*{h2l+h22\/^)Pp \h22?0l + \h2i+h22J^?Pp+0\ 



'p 

Similarly, when the CR user is active, the primary user's rate with given channel realizations 
can be computed as 

Rihiuhn)=I{Xp;Yp,Hn=hiuHi2 = hn) 



{l-ai)Pc\hnr + \y^hn + VaiPc-hi2r + oi^ 
^ {l-aMhn\^ + al " 

III. Parameter design for fading channels 

In this section we provide semi-analytical methods for finding the relaying ratios tti and 
precoding coefficients tti under fast and slow fading channels. 



A. Fast fading scenario 

With full CSIT, the CR transmitter must ensure the primary user's rate unchanged as in 
Viswanath's setting [4]. For fast fading channels, a meaningful definition of such rate is the 
ergodic capacity [7]. Therefore, we take the rate functions dT]) and dH) averaged over all channel 
realizations to form the following design criteria 

maximize E[R{H22,H2i)] (10) 
subject to £[i?(//ii,//i2)] >£[log(l+ '^"f^^ )]. (11) 

Since (fTTI) only depends on tti but not a2, we can solve (fTTI) for tti and then (fTOl) for CL2 
independently. The key to solving (fTTI) is to transform ^ into a ratio of quadratic form. 
By resorting to the moments of the quadratic form [11] we can find <X\ by solving a simple 
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logarithmic equation numerically. Thus we first rearrange the terms in ^ into the following 
matrix form. Let 

^p^(//ii,//i2r~C^(A^,E), (12) 
p = (iJiu^JiiV, and E = diag([ofi ofj])- ^1^° 

Q^(' ' ] andP-f Il_ 1. (13) 

\0 (l-ai)P, y \y^a^p aiPc 

Without loss of generality, we let both a|^^ and be 1. Then the primary user's rate formula 
becomes 

H^PH 

«(H,„H.)=.og(l + -^). (14) 

Note that the right hand side of (fTTI) is assumed to be a known constant Rg^g^^j^-. for the CR 
transmitter. Similarly, after some manipulations the CR user's rate becomes 

R{H22, H21) = log — rr-=V ^T^^^ +logO? , (15) 



where = [i/22,^^2i]^, cq = o?^ + |a2p/p, c = cq, and 
D = 



a^(a?^ + a2 v/^^T^)/';, |o| +a*yar^|2 



The optimization problem described in (flOl) and (fTTI) is not convex. Thus we propose methods 
described in Theorem [Hand Proposition [H to solve it sub-optimally with proofs given in Appendix 
lAl and El respectively. 

Theorem 1: The relaying ratio tti for the fast fading CR channel can be found sub-optimally 
by solving 

log(— -) - ■ ' ,2 = Krgodic^ (16) 



where 

A/e, =A/'SA/-ftr(ES), al=2iJ'SZS^i + tr{LS)\ ^i^, = A/'QA^ + tr(EQ), 
S = P + Q, /J and E are defined right after (fT2l) . and P and Q are defined in (fT3l) . respectively. 

After knowing ai, we can treat (fTOl) as an unconstrained optimization problem. Let B{a2) = 
(co(P + Q) - D) Hg + c, which is a convex function of a2. With the fact that the logarithm 
is concave, R{H22, Hii) is concave with respect to <X2- Therefore, from [?, p. 209] we know that 
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(flOl) is concave. Thus, we can set the derivative of E[R{H22, H21)] with respect to the complex 
conjugate of a2 to zero to find the extreme value 

where fn^ is the probability density function of H^., and the validity of the interchange of the 
integration and differentiation of the second equality is obvious. Note that solving the general 
moment of quadratic form [11] in (fTTl ) is intractable. Thus we resort to finding the suboptimal 
solution. 

Proposition 1: The precoding coefficient QL2 for the fast fading CR channel can be approxi- 
mated by 



(^l2/'2i + ,/¥^)(l-a,)P, 



^T^^, ■ 

B. Slow fading scenario 

For the quasi- static slow fading channels, the decoding error probability cannot be arbitrarily 
small since the transmitter does not know the reliable transmission rate with the limited channel 
knowledge. In this channel, the outage probability [7] for a certain target rate is more suitable 
than the Shannon capacity to measure the performance. 

For the CR user with a target rate R'"^, the outage probability which is defined as /'(i?(i/2 1,^22) < 
R^^) should be minimized. Let P^^^l and R^ be the primary user's achievable outage probability 
and outage capacity, respectively, in the absence of the CR transmitter. Since the CR system 
must ensure that the primary user's outage probability does not increase, we have the following 
constrained minimization problem for the outage probability of the CR user 

minimize Pf,^j 

subject to P{RiH2uH22) < R'^'') < Pfnft (19) 
P{R{Hn^H2i)<R')<P!;,,. (20) 

Assume that the CR transmitter knows the primary user's target outage probability Pf^^^J and target 
rate R^ . Similar to the ergodic case, tti and a2 can be solved from (|20l) and (fT9l ) sequentially. To 
reformulate the complex outage constraint (l20l) into an analytically solvable one, we resort to the 
concept of deterministic approximations of probability inequalities [8]. In [10] it was observed 
that if the mean of the rate is maximized, the deviation of the rate will also be increased, which 
may result in a worse outage capacity. As a result, if we want to achieve the outage probability 
constraint efficiently, we need to consider both the mean and deviation of the rate simultaneously. 
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In addition, to take the tightness of the approximation into consideration, we adopt the modified 
Cantelli's inequality [12]. The results are summarized as the following with proofs in Appendix 
O and respectively. The parameters /j^^ , P, and Q are defined in Theorem [TJ With 

Ai^ -^„P , (21) 

we have 

Theorem 2: The relaying ratio tti for the slow fading CR channel can be found by solving 
the following equation 

A^A, + y£^aA, = ^, (22) 

where /jai and are the mean and the standard deviation of Ai, respectively. The constant 
r can be 1 or 2/9 depending on the ^-factor. 

Since (|2TI) is a ratio of quadratic form, we can resort to [13] to efficiently compute the mean 
and the variance of Ai by the following proposition. Then (|22l) can be easily solved by a 
polynomial equation of tti. 

Proposition 2: The mean /jaj and variance in Theorem |2] are 

a / m s \ f t s 2m\ ^^^^ 
— 1 + — and— — + — , (23) 

respectively, where a = ;u"PAi + tr(EP), ^ = 4;u"PEP;U + 2tr(EP)2, t = 4A/"QEQAi + 2tr(EQ)2, and 
m 4 4;u"QEP;U + 2tr(QEPE) . 

Similarly, directly finding the optimal a2 from calculating the outage probability in (fT9l) is 
intractable. The main idea of the proposed method to find the semi-analytical solution of a2 
is described as follows. We first rearrange CR's rate formula i?(//2i,^22) into a weighted sum 
of non-central chi- square random variables which can be further approximated by a single chi- 
square random variable. Then we can semi-analytically find a2 minimizing the outage probability 
^oMf wi'^h the given outage capacity R'''^. This is summarized in the following theorem. Note 
that both v(a2) and w(a2) are simple rational functions of a2. 

Theorem 3: The precoding coefficient a2 for the slow fading CR channel can be found by 
solving 

/vv(a2) cd-l \ 

mm , , , (24) 

where d = {2^'^ — \ )/c5\^, T is the gamma function, y(.,.) is the incomplete gamma function, 
and 

, , ltr(EEEE)+2A/"EEEA/ , , , (tr(EE) +a^HEa/)^ 
^^^2) = 2 tr(EE)+^HE^ - "^^2^ = 'mEEEEHV^' ^'^^ 
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where E = (1 — cod)(P + Q) +dD, and c, cq are defined right after (fTSi) . 

The problem of solving (f24)) can be further approximated by the following corollary. 

Corollary 1: The precoding coefficient a2 for the slow fading CR channel can be approxi- 
mately found by solving 

w(a2) 

minfl-expf-^^^')')^, (26) 



V \ 2v(a2) 

where 

1, 0<w(a2) < 2 

(27) 

^ (r(l+w(a2)/2))^, w(a2)>2. 
Note that Corollary [T] does not need the computation of the incomplete gamma function in 

(|24|) . thus the search of a2 becomes simpler. 

IV. Asymptotic analysis 

Since a Rician fading channel becomes non-fading when its A'-factor approaches infinity, and 
the relaying ratio and precoding coefficient are well known for non-fading channels [4], we can 
verify the correctness of the proposed methods by making the /^-factor infinity. 

A. Asymptotic analysis of tti 

In the high ^-factor region, the channel gains are almost deterministic and knowing the mean 
values of the channels is almost the same as knowing the channel gains. With the assumption 
that + ojj = 1, 1 < J < 2, parameters and performances of slow and fast fading channels 
should converge to that of the non-fading channel with unity channel gain. Thus both Oej and 
approach zero. Besides, 

lim ;Ue, =£i, lim a/ez = ^2, 

K — >oo K — >oo 



where the definitions of ei and £2 are given right after (1431) . For the fast fading case, the left 
hand side of (fT6l ) becomes 

lim E[R{Hn,Hn)]=log{\±^)=log{\^)=log{l+A-')=\og{l +Pp), (28) 



where Ai is defined in (|2T]) . On the other hand, for the slow fading case (|22|) becomes 



lim ( + a/t^ - lOAi ) = -"A, = Ai =/'„ 1. (29) 



\ V out / 

From (|28l) and (|29l) we can find that tti converges in both fast and slow fading scenarios to the 
same value. From the definition of Ai we have 

Pp + aiPc + 2^aiPcPp 

(l-aOP. + l " ^' 
which results in the same tti as that in the non-fading CR case [4]. Thus we can conclude that 
the ai's found by the proposed methods in both fast and slow fading channels converge to that 
of the non-fading case when K 0°. 
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B. Asymptotic analysis of CL2 

Again the preceding coefficient (I2 should coincide with that of the non-fading case when 
K ^ 00. In the following we verify this statement. 

1) Fast fading: When K ^ 00, both ^22 and 1^21 approach 1, and (fTSi) becomes 

(1 + a/^)(1-«i)/'c 
iL-"^ - tk)P.,. = " ^ /f^)-?"-. (30) 

where \J^^ is due to the CR's relaying of the primary signal and a^^"^^ is the precoding 
coefficient of the non-fading case [4]. Thus it is verified that the proposed method also results 
in the MMSE linear assignment as the perfect CSIT case when K ^ 00, 

2) Slow fading: The asymptotic property of a2 is described in the following corollary with 
proof given in Appendix IB 

Corollary 2: When ^-factor — > 00, a2 derived in Theorem [3] converges to (1 + y'^^A)cx^^'^^. 

From (|30l ) and Corollary [2l we conclude that a2's for slow and fast fading channels converge 
to the same value in [4] when the A'-factor 00, 

V. Practical lattice-based linear-assignment Gel'fand-Pinsker coding 

All the previous results are based on the LA-GPC with theoretical, unstructured Gaussian 
random codebooks. In the following we introduce a practical lattice-based coding to implement 
the LA-GPC for the CR user. With the carefully designed filters specified in Theorem H the 
proposed lattice-based LA-GPC is a non-trivial extension of the well-known lattice-based DPC 
[14]. First, the side-information filter is selected according to the aforementioned precoding 
coefficient a2 instead of directly using the MMSE filter ttc described after Q for the interference- 
free channel as in Section lTl-Bl of [14]. Secondly, with only the statistics of the CSIT, the receiver 
filter selection is more involved than simply choosing it as ttf in DPC [14] [15]. Specifically, 
this filter must be different from the side-information filter, and this is contrary to the common 
practice in the lattice-based DPC [14] [15]. Finally, we adopt the finite dimensional self-similar 
lattices to implement the lattice coding structure which is much more feasible than the lattice- 
based DPC in [14] using a very long codeword length. The proposed coding works very well 
with a reasonable codeword length (and decoding latency), and our simulation in Section IVI-CI 
verifies this. 

To illustrate the lattice-based LA-GPC, we focus on the following channel corresponding to 
^ with discrete time index t as 

Y{t) = H,{t)XL{t) +Hs{t)S{t) +Z(0, (31) 
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where I <t <T, T is the codeword length. To emphasize the difference between lattice coding 
and random Gaussian codebooks used in the previous sections, we use Xi(t) instead of X{t) to 
represent the lattice-coded signal. As in Section III-C[ we first focus on the rate Q with given 
channel realizations Hx{t) = hx{t) and Hs{t) = hs{t). We also rewrite (|3TI) in an equivalent real 
super channel to present the lattice coding more easily. By concatenating all T symbols, the 
channel becomes 

y =HxX^ + Hs5 + Z, (32) 

where = (Xj j, . . . ,XJt)^ and X^, = {KQ{XL{t)} ,\m{XL{t)]y . The non-causally known side- 
information vector at the transmitter S_ and the noise vector Z are obtained similarly from S{t) 
and Z{t), respectively, as from Xt{t). The covariance matrices of 5 and Z are denoted by 
and E^, respectively. The 2T x 2T block-diagonal real channel matrix Hx is diagi^x^ ■ ■ -Hx"^]), 
where the tih. diagonal term is 

\R^{h,{t)} -\m{h,{t)}' 
' [lm{h,{t)] Re{/i,(0} J' 
and Hs is formed from Hs{t) in the same way as Hx from Hx{t). The channel input power 
constraint is (1 —ai)Pc/2 because the CR user transmits its own signal with power Pc/2 per 
real dimension, and relays the primary user's signal with power <X\Pc/2. 

We will first give a brief review of the lattice codebook, then introduce our proposed lattice 
coding. A 2r-dimension real lattice A is defined as A = {Gb : b G 1?^}, where G is the 2T x 2T 
generator matrix of A. The Voronoi region 'V is the set of points g G R^^ which are closest to 
in Euclidean distance than to other lattice points X G A. Every g G R^^ can be uniquely written 
as g = X-|-ng, where X G A and G i^. With quantizer input g, the lattice quantizer associated 
with 'P is defined as 2(g) = X, if g G A, + The modulo-A operation associated with is then 

gmodA = g-2(g). (33) 

We define the nested-lattice codes with length 2T as 

Definition 1: Let Ac be a lattice and A^ be a sublattice of it, that is, A^ C A^. The codeword 
set of the nested lattice code is Cc = {A^ mod A^} = {Affl The code rate of this nested 
lattice code is Rc = ^log | l^^^l |/| |'^'c| |, where 'Uc and i/"^ are the fundamental Voronoi regions of 
Ac and A^, respectively. And for a bounded Jordan-measurable region n/ C M?'^ , ||t^|| denotes 
the volume of . 

For conciseness of the paper, we focus on the slow fading channels where the channels remain 
constant within a codeword length. The proposed scheme can be easily modified for the fast 
fading channels [16]. The proposed coding works as follows 

Encoder: The encoder selects a codeword Cc G Cc according to the message index and sends 

X^ = V(l-ai)Pc((Cc-FsS-D) modA^) , (34) 
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l7 



(35) 



where the IT x IT side-information filter Fs is formed from <X2 in Theorem |3] as 

Re{a2} -Im{a2} 
Im{a2} Re{a2} 

and ® denotes the Kronecker product. The dither signal D, uniformly distributed in "Vq and 
independent of the channel, is known to both the transmitter and the receiver. This dither plays 
a critical role in making independent of Cc and 5 [16]. The second moment P{i'q) [15] of 
Aq is set 1/2 to satisfy the power constraint (1 —a\)Pc/2. 
Decoder: The decoder performs the following operation on Y_ 

£ = L(Fry+D). (36) 

After some manipulations [16], £ in (l36l) can be rewritten as 

F = L«)+e, (37) 

where G A^ + Cc C Ac and 

e^ (FrH-l2r)^+(FrHs-F,)S + FrZ, (38) 

where H = a/(1 — ai)/'cHv. The receiver filter Fr is obtained according to the MMSE filter for 
estimating the auxiliary random variable U in ^ from 7 in ([21) as [16] [17] 

Fr = (^H^ + WE,hJ)(^HHT + HsE,HsT + E,)-i, (39) 

where E^ is the covariance matrix of 5. The whitening matrix filter L for e must satisfy 

l7l. = Y.-^\ (40) 

where E^ is the covariance matrix of e when T — oo as 

He = i(FrH-l2r)(FrH-l2r)^ + (FrH,-Fs)E,(FrHs-F,)T + FrE,F/. (41) 

Then we can use the generalized minimum Euclidean distance lattice decoder to decode Cc. 
First the decoder finds 

b = argmin|y-LGcb|2, (42) 

where Gc is the generator matrix of the channel coding lattice Ac. After that the decoded 
codeword is Cc = [Gcb] mod A^. 

With "good nested lattices" defined in [15], we have the following result. The proof is omitted 
and can be found in [16]. 

Theorem 4: Let the side-information filter Fg, receiver filter F^, and whitening filter L be 
selected as (|35l) . (|39l) . and (|40l) . respectively. Based on sequences of good nested lattices, the 
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coding specified in (I34l)-(|42|) with ([3]) is able to achieve the linear-assignment rate i?(/i225^2i) = 
Y^og l^^l^^ under power constraint (1 —ai)Pc when r — > oo where is defined in (|4TI) . 

It can be proved that with full CSIT, the proposed coding reduces to the DPC in [14] [15]. 
And according to our simulation, the DPC suffers severely when it is naively applied to channels 
without perfect CSIT. 

VI. Simulation results 

In this section we demonstrate the performance of the proposed LA-GPC based CR system 
under both fast and slow fading channels. We also show the performance of the proposed practical 
lattice encoder/decoder which implements the LA-GPC. For simplicity, the four channels in Fig. 
\T\ are assumed independent Rician with the same ^-factors. 

A. Fast fading 

To verify the goodness of the proposed approximation, we first compare it with brute-force 
full search in terms of the primary user's rate. Signal and noise powers are set as Pc = Pp = 10 
and o|^^ = = 1. This comparison is shown in Fig. [2l The unit of the vertical axis is bits per 
channel use, bpcu. We can see that the proposed approximation performs well for all ^-factors. 
That is to say, the primary user's rate is approximately the same as that with tti obtained by 
brute-force search. More specifically, the proposed method may over-design ai slightly at the 
small ^-factor region such that the primary user's rate becomes larger than the target. At the 
large ^-factor region, ai obtained by the proposed method almost coincides with that obtained 
by full search. The CR user's ergodic capacities using tti and tti from the proposed method is 
shown in Fig. [3l It can be observed that the over-designed ai slightly increases the primary's 
rate at the cost of the CR user's rate. However, since ai is only slightly over-designed, the 
resulting change of the CR user's rate is small. In addition, we also compare cases of treating 
interference as noise and with full CSIT. The former neglects any useful information to cancel 
the interference, and incurs a huge rate loss for all ^-factors considered. The ergodic capacities 
for the partial and full CSIT cases converge when the ^-factor is large. This phenomenon has 
been verified in Section |IVl The rate loss at the low ^-factor region is due to two factors. One 
is the imperfect CSIT. The other is that the linear assignment strategy may not be optimal for 
the partial CSIT case. Until now, the optimal strategy for the partial CSIT case is still an open 
problem. Note that all the LA-GPC curves can also be obtained from the lattice-code achievable 
rate l/TloglE^I^^ with filters specified in Theorem |4] and T ^ oo. The naive DPC scheme 
assumes that the channels always take values at their means, and uses these mean values to 
design the DPC filter. That is, the Fs filter is naively chosen according to ttc in Section III-BI 
with Hj^ replaced by E[H^, and Fr = FsHx^^ It can be observed that the naive DPC scheme 
has a significant rate loss compared to the proposed method in a wide range of K-factors. As 
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expected, the ergodic capacities of all cases except treating interference as noise converge to the 
same value when the ^-factor is large. 

B. Slow fading 

We first demonstrate the effectiveness of the approximation for ai in terms of the primary 
user's outage capacity. A good tti will make P{R{H\i,H2i) <R^) as close to P^^^^ as possible. 
We then compare the values of a2 obtained from (|24|) and full search. Finally, we illustrate the 
CR user's outage probabilities with a2's from these two methods. Signal and noise powers are 
set as Pc = Pp = 10 and o^^ = = 1. Since different channel conditions can support different 
outage capacity and outage probability pairs, it is more reasonable to consider the supportable 
pairs than using the same pair for different ^-factors. Thus, in the following, we set primary 
user's outage capacity (in bpcu) and probability pairs as (1, 0.1), (2, 0.1), (2, 0.01), and (2, 
0.01) for ^-factors of 0, 5, 10, and 15 dB, respectively. Moreover, the CR user's target rates 
are set as 0.2, 0.5, 1, and 1.5 bpcu, respectively. Fig. |4] checks whether the primary user's 
performance constraint (|20l) is satisfied. Compared to the baseline case without the presence of 
the CR transmitter, we find that the approximation of tti is very tight in the low ^-factor region, 
and is good for medium and high ^-factors. The tight approximation is a result of the ^-factor 
dependent r. When the ^-factor is or 5 dB, it can be found that 5 > and the distribution of 
^1 /Ai is asymmetric. Thus we let r = 1. On the other hand, when the ^-factor is 10 or 15 dB, 
/Ai is unimodal and symmetric simultaneously. In addition, 6 > 2/v^. Thus we set r = 2/9 to 
further tighten the original Cantelli's inequality. The slightly over-designed ai ensures that the 
primary user's outage probability constraint is met at the cost of the CR user's rate. 

The CR user's outage probabilities with tti and a2 both obtained from approximations and 
full searches are given in Fig. \5\ Recall that ai is slightly over-designed at medium and high 
A'-factors. This reduces the CR's power for its own signal, and increases the outage probability. 
On the other hand, when the ^-factor is small, the worse channel condition dominates the 
performance. In addition, we can find that the approximation from Alzer's bound [?] of the 
incomplete gamma function is good for the considered ^-factors. Naive DPC again incurs a 
large rate loss in the low to medium ^-factor region. It can be seen that if the interference is 
treated as noise, the outage probability is much worse than those of the proposed precoding 
methods. Similar to the fast fading case, there is a gap between the outage probabilities of the 
full and partial CSIT cases. Since the optimal precoding strategy is still an open problem, the 
minimum possible gap is unknown. The non-decreasing outage probability curves of treating 
interference as noise and full CSIT are due to the facts that the considered target rates are not 
identical for different ^-factors. 
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C. Lattice implementation with finite codelength 

We investigate the error performance of lattice precoding with a reasonable codeword length 
(and decoding latency) at large enough SNR. A pair of self-similar nested lattices is used, with 
the Gosset lattice Es [18], which has the densest packing in 8-dimension, as the fine lattice 
Ac. The reason to use the Gosset lattice Eg is that in moderate dimensions, it is well-known 
that the best lattices in terms of coding gain are also good quantizers. That is, the lattice can 
also be good A^. The coarse lattice A^ is generated by A^ = QAc where the coding rate is 
given hy R = 21og2. Recall that transmitted signals is assumed to be Gaussian in previous 
derivations. The Gaussianity of the signals at the proposed lattice-based CR transmitter output 
with a reasonably short codeword length T = 4 is shown in Fig. [6l It is clear that the transmitter 
output distribution is almost Gaussian. 

In the simulation of codeword error probabilities we consider two ^-factors: OdB and lOdB 
and two code rates: 2 and 4 bits per channel use. The noise variance is normalized to 1 and the 
primary user's power constraint is assumed to be Pp = 100. The LA-GPC curves are obtained 
from the lattice-code achievable rate with filters specified in Theorem |4l A Fano sequential- 
decoding based lattice decoder [19] is used to solve (|42|) . We compare the lattice coding results 
with the theoretical results with random Gaussian codebooks. The basis for the comparison is that 
the codeword error probability is approximately the same as the outage probability. The results 
are shown in Fig. |7] and Fig. [S] For comparison, we also show the results without interference 
from the primary user and treating interference as Gaussian noise. When the interference is 
treated as noise, the codeword error rate is 1 regardless of the SNR and ^-factor ranges. It can 
be seen that in the considered A'-factor range, the performance of the proposed nested lattice 
decoding is approximately the same as that derived theoretically. This verifies the goodness 
of the used Gosset lattice code. We also observe that in Fig. |7] with ^-factor = OdB, lattice 
precoding with partial CSIT performs slightly better than the theoretical outage result. A similar 
phenomenon was also reported in [20]. 

VII. CONCLUSION 

In this paper we considered the cognitive radio channel with partial CSIT. Using the linear- 
assignment Gel'fand-Pinsker coding, we proposed semi-analytical methods for finding the re- 
laying ratios and the precoding coefficients for both fast and slow fading channels. Asymptotic 
analysis showed that the relaying ratios and the precoding coefficients obtained by the proposed 
methods converge to those with full CSIT when the ^-factor approaches infinity. Simulation 
results showed that the proposed semi-analytical parameter design methods perform almost 
as well as exhaustive search. We also used nested-lattice coding and decoding to realize the 
precoding in practice. Simulation results showed that the nested-lattice coding scheme can 
effectively reach the achievable rate of the CR channel with partial CSIT at high SNR. 



16 



Appendix 

A. Proof of Theorem [7] 

After applying the expectation operation to (fT4l) we have 

E[R{Hn,Hn)] =E[\og{l +^^^8^^)] -£[log(l +^?Q^p)] 

^£[log(l+£i)]-£[log(l+£2)], (43) 

where S = P + Q, Ei = H^SH_p, and £2 = KpQKp- To guarantee that (fTTI) is valid, we expand 
/(Ei) = log(l +ei) and ^(£2) = log(l +£2) by the Taylor series with different orders, respectively. 
The reason is as the following. We first expand log(l +£1) by the k-th order Taylor series respect 
to the mean of £1, /jgj, as 

log(l+£i)=log(l+^e,)-log.Et2^^^^^ + oi^A(£l)+oi, (44) 

where fk{£i) is the truncated k-th order Taylor expansion and oi is the Lagrange remainder. 
After applying the expectation operator to (|44l) and truncating the remainder o we have 

Em = iog(i -iog.EL2 ^'5^nV%^"^ - (45) 

where E[fj^{ei)] is further simplified as E[fj^]. Since the Taylor expansion of the logarithmic 
function / is an alternating series, and £[(£1 -A'ei)"+^] > 0, Vn G IN, from [11], we have the 
following relations: E[fi] > ^[/g] > E[f5] > ■ ■ ■ > E[f] and E[f2] < Eif^] < E[fe] < ■ ■ ■ < E[f]. 
Similar properties can be found for ^(£2)- The truncated first order Taylor expansion of £'[log(l + 
£2 )] can be derived in the same way. From the above we know that E[f]—E [g] > E [f2i] — E [gij- 1 ] ■ 
Thus solving the following equality 

E[f2i] -E[g2j-l] =Rergodic 

makes the constraint (fTTI) valid where i,j E IN. Here we choose i = j = 1 for simplicity. After 
some manipulations we have (fT6l ), where /Jej, Og^, and /jg^, the mean, variances of £1, and the 
mean of £2, can be found according to [11]. 



B. Proof of Proposition [7] 

With the fact that co(P + Q) — D is positive semi-definite and c > 0, B{a2) is positive. Then 
we can generalize the result in [21], which is aimed for the ratio of quadratic forms, to the ratio 



of general quadratic formJll in (fT7]) . Therefore, we have 



'B\a2) 
_B{a2) 



E[B'{a2)] 
E[Bia2)] 



0. 



"General" means that there is not only the term of the quadratic form, but also terms of lower orders. 
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Since B{a2) is a second order polynomial of a2, to proceed, we only need to consider the 
numerator, that is, 

E[B' {a2)]=0. 

With the fact that O22 + |A'22p = 1, the resulting precoding coefficient with the statistics of CSIT 
is as (fT8l). 



C. Proof of Theorem |2] 

From (O, (l2Ql) and dZB, the optimal tti will make 

/'(Ai> ^^)=/'f„,. (46) 

To make this problem analytically tractable, we relax the equality constraint by inequality to 
find the sub-optimal solution. With the non-negativity of the transmission rate, the Cantelli's 
inequality [22] is a good candidate for finding the solution. Considering additional properties of 
the distribution //Xj, such as the modality and symmetry, the Cantelli's inequality can be further 
tightened [12]. Thus we propose the use of the following modified Cantelli's inequality to relax 
(l46l) as 

P(Ai>a/a,+5oaJ< (47) 

According to different resulting from different channel conditions, r in (|47]) can be catego- 
rized as [12] 

' 1, if5>0, 
§5 if S > and/Ai is unimodal and symmetric. 

Comparing the left hand side of the inequalities in (|46|) and (|47]) results in 

1 

2^"~ 1 

Also by comparing the right hand side of the inequalities in (|46l ) and (1471 ) we have 5 = 
a/tT^]^^^^. Substituting 5 into (|49l ). we have (|22|) . To determine r, we can simply substitute 
r=lorr = 2/9to see whether the resulting 6 meets the constraints (|48l) . The constraints 
of modality and (or) symmetry which may be satisfied under different channel conditions are 
discussed in Section |VIl Thus we can solve the outage probability constraint by finding <X\ 
which satisfies the equality in (|22l) . We now discuss the validity of (|20l ). As for the effect of the 
modified Cantelli's inequality, note that the modified Cantelli's upper bound is larger than the 
real probability value [12]. Since we equate the upper bound to the target outage probability to 
solve tti, the resulting ai will be larger than the optimal one from (|46l) . 



//A1+6OA1 =;;^P — 7. (49) 
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D. Proof of Theorem 13 

First, substituting (fT5l) into (fT9l) . and after some manipulations we have the following opti- 
mization problem 

min P{z = H^EH^- = ^2 £2^2 < - 1 ) , (50) 

where E2 = E^/^EE^/^, and H2 ~ C9i {L^^/^iu.h). The transformation from E to E2 is to make 
H2 have an identity covariance matrix for the convenience of the following analysis. After 
substituting the eigen decomposition E2 = v^Av into z, we have 

z = HfAH^ ~ Xix2(2, 2|a/3i I') + W(2, 2|a/32 P) , 

where H_2, ~ C5V; (/J3,l2), /J3 = [1^31, — vE^^/^/j, and Xi and X2 are the eigenvalues of E2. 
The notation X^iP:^) denotes the noncentral chi-square distribution with degree of freedom p 
and non-centrality q. Note that the objective function (l50l) is a function of a2. That is, Xi, X2, A'si, 
and ;U32 are controlled by a2. 

To calculate the outage probability (|50l) . we resort to approximating z by a scaled single central 
chi-square random variable [23] as 

z~^a2(2,2|A^3i|')+W(2,2|A/32|')^v%2(w). (51) 

Recall that the chi-square distribution is a special case of the gamma distribution. Thus w does 
not have to be an integer. The parameters v and w should be chosen such that both sides of (|5TI) 
have the same first and second moments [23] as 

vw = Xi(l + IA/31P) +?t2(l + |a'32|^), (52) 
2v^w = (1 -f2|A/3iP) +Xi(l +2|A/32p). 

With the fact 

;ii+?i2 = tr(E2)=tr(EE), (53) 

h + = pfAiJ3 = a/"E^v«AvE^;/ = iJ^E^, (54) 

we can find the value of vw. Similarly, from 

Xl + Xl = tr(E2E2) = tr(EEEE) (55) 

Allwi P + ^IM^ = A/"E^E2E2E^A' = a/"EEEa/. 

we can find the value of 2v^w. Then v and w can be expressed as (|25l) . After that, solving (l50l) 
with (I5TI) is equivalent to solving 



cd-l 



mm 



which can be further represented as (|24l) . 



19 



where 



E. Proof of Corollary \T} 

From [?] we know that the incomplete gamma function can be bounded as following 

r(a) 

; ^' 0<a<l 
s = < -I (56) 

^ (r(l+a))-, a>l. 

Let a = ^^^^^ and x = we then have (|26l) . 

Proof of Corollary |2] 

As the ^-factor— > 0°, the distribution of z defined in (|50l) becomes a delta function with a 
nonzero value at its mean fj^. In the following we need to find a2 maximizing the mean value 
of z which can be approximated by vw, as shown in (ISTI) . Since c in (l50l) is a function of a2, 
we must take it into consideration. The problem becomes to maximize ^ cJ. As both /J21 and 
A^22 ^ 1. and and ^ at very high ^-factor, from (|52|) . (|53l) . and (|54l) . we have 

limpz = vw = iLi^EiLi + tr(EE) ~ en + ei2 + ^21 + £22, 

K — > 

where en, ei2, e2i, and e22 are the four entries of E. Since i^z — cd is a quadratic formula of a2, 
we can take derivative to find the maximum as following 

1 a 



d dal 



ilj,-cd) = -{{l-ai)Pc+l) Ppa2 + {Pp + ^aiPcPp) ( 1 - ai )/',. = 0. 



Then (I2 can be found as 



l + J^){l-ai)Pc 



which is the same as that of the non-fading case. 
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Fig. 1. The model of the cognitive radio channel. 
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Fig. 2. Comparison of the primary user's ergodic capacities under full search and the proposed method with Pc = Pp = 10 and 
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Comparison of the CR user's ergodic capacities under full search and the proposed method with Pc = Pp = 10 and 
. = 1. 
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Fig. 4. Comparison of the primary user's outage capacities obtained by the proposed method and the target rates with 
p^=Pp = 10 and = = 1. 
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Fig. 5. Comparison of the outage probabilities with a2 from full search and the approximation with Pc — Pp = 10 and 
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Fig. 6. The Gaussianity of the CR transmitted signal using the nested lattice code. 
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Fig. 7. Comparison of tlie outage probabilities with rate= 2 bpcu. 
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Comparison of the outage probabilities with rate= 4 bpcu. 



